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1 Introduction 

The work reported herein represents a continuation of work begun at MSFC during the 1991 NASA/ASEE 
Summer Faculty Fellowship Program [5, 7]. During this period the author proposed and studied a paradigm 
for the analysis and synthesis of Integrated Health Monitoring and Control Systems (IHMCS) for rocket 
engines. This work was motivated by NASA’s desire to develop advanced propulsion systems which could 
operate with increased reliability at decreased cost. The Space Transportation Main Engine (STME) whose 
function it is to provide primary thrust for the National Launch System (NLS) family of vehicles is a primary 
example of this trend [2, 9]. 

The need for an integrated approach to health monitoring and controls was established by the author 
who has identified significant interaction effects which exist between control and health monitoring functions 
[5, 7]. The nature of these interaction effects is such that unless they are taken into account during the 
design phase, significant performance degradation may occur in overall the system operation. 

The specific approach proposed allows the designer to embed the IHMCS into a general system archi- 
tecture wherein the wide array system analysis and design tools can be brought to bear. Within the this 
approach off-nominal conditions are modeled as indicated in Figure I. Here p denotes some nominal com- 
ponent or subsystem within the rocket engine. In the first case (Figure 1(a)), off-nominal conditions are 
represented by exogenous signals / injected at either the component input or output. When off-nominal 
conditions of this sort are incorporated into the overall engine model, the IHMCS analysis and design tasks 
can be reduced to problems in tracking and disturbance rejection [5, 7]. In the second case (Figure 1(b)), 
off-nominal conditions are represented by exogenous component dynamics A which may augment the nomi- 
nal component dynamics. When off-nominal conditions of this sort are incorporated into the overall engine 
model, the IHMCS analysis and design tasks can be reduced to problems in uncertainty accommodation and 
robustness [5. 7]. 

In either case, it is clear that the approach discussed above is model-based. Thus, in order to apply 
these results it necessary to have access dynamical models for both the nominal rocket engine system and 
all off-nominal conditions of interest. Alternatively, and perhaps more favorably, would be to have a method 
for obtaining such models in a systematic manner. It is the task of this report to present just such a method. 

This methodology is based on the application and manipulation of the fundamental laws of conservation 
in order to derive dynamical models for thermo-fluid systems such as those contained in chemical propulsion 
systems. It has a number of significant features which make it well suited for use addressing the problems 
encountered in IHMCS design: First, it allows the assembly of dynamical nominal and off-nominal engine 
models of low order by indicating the significant dynamics. Second, it allows us to distinguish between those 
off-nominal conditions which are likely to be modeled as signal type (Figure 1(a)), and those likely to be 
modeled as uncertainty type (Figure 1(b)). Third, it allows for the easy incorporation of various sensor and 
actuator types. Finally, it provides models in a format suitable for direct computer simulation. 

Due to space constraints, the development here is necessarily brief and the interested reader is referred 
to [6] for a more detailed treatment. 


2 Development of Thermo-fluid Modeling Principles 

Chemical propulsion systems such as the STME are basically thermo-fluid systems whose dynamics are 
governed by the laws of conservation of mass, momentum, and energy. Discussion of these basic laws of 
physics can readily be found in the literature (e.g., [3, 8]). Here we will apply these laws in order to develop 
a set of principles for model development by studying the dynamical behavior of the fluid in a generic engine 
component represented by the variable area control volume given in Figure 2. The results given here are 
motivated by treatments given in [8, 4]. 

The basic assumptions which are made for the purposes of the development below are given as follows: 
(i) Fluid is allowed to cross the boundary only at the inlet (x = xi) and the outlet {x = xo). (ii) We assume 
that the fluid flows in the component are quasi-one-dimensional, i.e., flow properties are axisym metric and 
uniform in any plane normal to the direction of flow, (iii) Where the working fluid (i.e.. the propellants, 
combustion products, etc.) are in gaseous state, they behave as ideal gases. Specifically, p = pRT , e = c v T< 
and h = c p T. (iv) Body forces are negligible, (v) The flow is inviscid. 
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Given these assumptions we can now develop dynamical expressions characterizing the behavior of the 
fluid in the generic engine component of Figure 2. 1 The behavior of this fluid can be completely characterized 
in terms of 4 state variables: V(x,t) - flow velocity, p(xj) - pressure, T(x,t) - temperature, and p(x,t) - 
density. These 4 variables are constrained to obey the laws of conservation of mass, momentum, and energy 
which can be expressed for the control volume of Figure 2 as follows: 

Conservation of Mass: 



Here the terms, from left- to- right, can be interpreted as follows: the rate of increase of mass in V, and the 
mass flow out of S. 

Conservation of Momentum; 

1////™ + ~/jG* (2) 

Here the terms, from left-to-right, can be interpreted as follows: the rate of increase of momentum in V, the 
momentum flow across S , the total body force in V, and the total pressure force on S . 

Conservation of Energy: 


° +W **-II s pVds = &JJl f {' + T) dv + JJ/{ e + T) vds W 

Here the terms, from left-to-right, can be interpreted as follows: the rate of heat added across 5, the rate of 
shaft work done in V , the rate of work done on S by pressure forces, the rate of change of energy in V. and 
the rate of flow of energy across 5. 

Compressible Flow: We begin with the case where the fluid flow is considered compressible. The as- 
sumption of compressibility is necessary in the rear stages of a rocket engine system where the propel- 
lants/combustion products are in gas phase. Under these conditions equation 1 can be simplified to yield: 


dp _ A 
dt ~ V 




where A denotes the mean area between inlet, x 1t and outlet, 

It is clear from this expression that if the ratio A/V is large, then the dynamics corresponding to the 
mass equation for the component tend to equilibrium quickly. In such cases, the dynamic mass equation can 
be replaced by the algebraic equation Pi^i^i = p,V 2 A 2 for that component, and the mass equation does 
not contribute to the system’s dynamic order. 

In a similar manner, equation 2 can be simplified yielding the expression 


djpV) 

dt 


1 


*2 ~ Xi 



T / *> A 1 

- 2 — + Pl — 



From this expression it is clear that if l/(x 2 - r,) is large, then the dynamics corresponding to the 
momentum equation for the component tend to equilibrium quickly In such cases, the dynamic momentum 
equation can be replaced by the algebraic equation Pl V^Ai + Pl .4i = p 2 V 2 2 A 2 + p 2 A-, for that component 
and the momentum equation does not contribute to the system's dynamic order. 

Finally, equation 3 can be simplified yielding the expression 


d_ 

dt 




{P\d\Vic P Ti - p 2 A 2 V 2 c p T 2 + W shaft + Q ) . 


From this expression it is clear that if V is small, then the dynamics corresponding to the energy equation 
lor the component tend to equilibrium quickly. In such cases, the dynamic energy equation can be replaced 
by the algebraic equation Pl A l V 1 c p T 1 + W shaft + Q = p 2 A 2 V 2 c p T 2 for that component, and the energy 
equation does not contribute to the system's dynamic order. 

Table 1 provides a catalog of the symbols used here. 
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Incompressible Flow Dynamics: Next, we consider the case where the fluid flow is incompressible. The 
assumption of incompressibility is especially reasonable in the front stages of a rocket engine system where 
the propellants are still in liquid phase. In this case p{x,t) is constant and so not a dynamic state variable, 
and the energy equation is not needed. The equations for conservation of mass and momentum given above 
can be simplified to yield 

Summary: Based on the analysis given here, the relative sizes for 


A 1 l A 
V ’ X 2 - Xi * V p(x 2 - Xi) 

for a given engine component can be used to decide which dynamics are required to model that component, 
and which dynamics can be replaced by algebraic relationships. 


3 Application to the STME 

The results outlined above were applied to develop preliminary models of the STME for the purposes of 
studying the health monitoring and control functions. 

Figure 4 provides a schematic diagram of the engine [2]. As it indicates the STME will employ a gas 
generator cycle with liquid hydrogen and oxygen propellants. Figure 3 provides a preliminary indication of 
the physical layout of the STME roughly to scale thereby providing of the relative dimensions of the various 
components [2]. 

Based on the information contained in Figure 4 the STME is modelled using the concepts from Section 2 
by first breaking it up into four major assemblies: Main Combustion Chamber, Fuel Turbo- pump, Oxidizer 
Turbo-pump, and Gas Generator. Next, the relative dimension information given in Figure 3 is used to 
characterize the behavior of each subassembly using the appropriate combination of dynamic and algebraic 
equations. This resulted in a 17th order dynamical engine model. Finally, the overall model was encoded 
into MARSYAS [1] and used for simulation studies. 
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Figure 1: (a) Signal type and (b) Uncertainty type off- nominal conditions. 



Figure 2: Generic engine component control volume. 

Symbol Variable 
p pressure 

e internal energy 

h enthalpy 

p density 

T temperature 

V flow velocity 

-4,5 area, surface area 

V volume 

Q heat 


Table 1: General variable definitions. 
Lti?- lex 



Figure 4: STME flow schematic. 3: S ™ E phyS1Cal lay ° Ut - 
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